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Abstract 


A  program  is  presented  for  the  generation  of  tables  of  solar  depression 
versus  time  on  a  particular  day  or  period  of  days  for  any  location.  The  input  data 
required  are  two  values  for  each  day  and  the  value  of  AT  for  the  year  as  taken  from 
"The  American  Ephemeris  and  Nautical  Almanac, 11  the  longitude  and  latitude  of  the 
location  of  interest,  the  maximum  and  minimum  zenith  angle  desired.  The  program 
is  suitable  for  or  easily  adapted  to  any  first  generation  or  newer  computer  capable 
of  being  programmed  in  Fortran. 
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I.  Introduction 

The  tabulated  values  of  solar  depression  versus  time  currently  available  are 
of  questionable  value  when  used  as  a  common  denominator  in  the  reduction  of  data 
gathered  during  atmospheric  observations.  In  most  cases,  the  discrete  intervals  of 
depression  angle,  day  number  and  latitude  are  too  widely  spaced  to  allow  for  an 
accurate  interpolation  between  given  values.  The  accurate  calculation  of  the  time 
associated  with  a  given  depression  angle,  at  a  specific  location,  on  a  particular  day 
or  period  of  days  has,  up  to  now,  been  a  laborious  task.  The  program  presented  in 
this  report  has  been  developed  to  allow  the  generation  of  useful  tables  for  any  loca¬ 
tion.  The  input  data  required  are  two  values  for  each  day  and  the  value  of  AT  for  the 
year  as  taken  from  "The  American  Ephemeris  and  Nautical  Almanac"  which  is  pub¬ 
lished  in  advance  for  any  year,  the  longitude  and  latitude  of  the  location  of  interest, 
the  maximum  and  minimum  zenith  angles  of  interest  and  the  increment  desired 
between  the  given  zenith  angles. 

II.  Program  Explanation 

The  basic  method  employed  for  computation  is  that  illustrated  in  section  13.  D 
of  the  "Explanatory  Supplement  to  the  Astronomical  Ephemeris  and  Nautical  Almanac." 
As  can  be  seen  upon  checking  the  sample  calculations,  reproduced  in  Appendix  E, 
the  accuracy  of  this  method  is  greatly  influenced  by  the  selection  of  the  time  used  for 
the  interpolation  of  the  tabulated  values  of  apparent  declination  and  ephemeris  transit. 
These  tabulated  values  are  taken  from  the  Sun  data  for  "0"  hours  ephemeris  time  as 
given  in  "The  American  Ephemeris  and  Nautical  Almanac, "  hereafter  abbreviated 
"AE  and  NA'.' 

In  order  to  achieve  maximum  accuracy,  a  reiterative  method  of  calculation  is 
used  based  upon  the  Newton  Forward  method  of  interpolation.  An  approximate  time, 
used  as  an  interpolation  factor  is  inserted  and  with  this,  a  calculated  time  for  the 
zenith  angle  of  interest  is  generated.  This  calculated  time  is  then  compared  with 
the  approximate  time  inserted,  and  unless  the  difference  between  the  two  is  less 
than  a  specified  tolerance,  (0.0001  hours),  the  approximate  time  is  replaced  by  the 
calculated  time  and  the  procedure  repeated  until  the  tolerance  condition  is  met.  Note 
that  this  leaves  us  with  a  value  of  ephemeris  time  for  the  angle  of  interest.  The  value 
of  AT  is  now  subtracted  from  the  final  calculated  time  to  yield  Greenwich  mean  time 
at  the  Greenwich  meridian.  To  this  value  is  added  the  value  of  West  Longitude,  con¬ 
verted  to  hours  to  finally  yield  the  local  Greenwich  mean  associated  with  the  particu¬ 
lar  zenith  angle.  The  program  next  increments  to  the  next  desired  solar  angle,  and 
the  procedure  repeated  until  all  desired  angles  for  a  given  day  are  computed.  Fol¬ 
lowing  this,  the  data  for  the  next  day  are  read  in,  and  the  calculations  are  repeated. 
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The  zenith  angle  used  to  calculate  sunrise  and  sunset  is  90°,  zero  degrees  hori¬ 
zon  angle,  rather  than  90°  50'  as  used  in.  the  Explanatory  Supplement  to  the"AE  and 
NA."  This  angle  was  chosen  because  the  primary  interest  was  in  the  center  point  of 
the  sun's  disk  rather  than  the  upper  limb;  therefore,  the  semi-diameter  of  16'  was 
neglected.  The  horizontal  refraction  of  34'  was  neglected  because  our  interest  was 
in  the  effect  of  the  sun  upon  the  upper  atmosphere. 

The  horizon  angles  that  are  tabulated  in  the  output  are  those  for  points  at  sea 
level.  Since  our  interest  was  in  zenith  atmospheric  measurements  only,  this  suited 
our  purpose.  If,  however,  the  altitude  of  the  observing  location  is  to  be  taken  into 
consideration,  the  zenith  maximum  and  minimum  angles  may  be  increased  by  a  fac¬ 
tor  of  1. 17h^“,  where  h  is  measured  in  feet. 

The  basic  formula  used  for  computation  is: 

cos(h)  =  -tan(6)  tan(cp)  +  cos(z)  cos(6)  sec(cp) 

where  op  is  the  latitude  of  interest, 

6  is  the  apparent  declination  of  the  sun  at  the  time  of  interest, 
z  Is  the  desired  angle  of  the  sun  as  measured  from  zenith, 
and  h  Is  the  local  hour  angle  of  the  sun. 

In  temperate  latitudes,  no  check  normally  need  be  made  of  the  value  of  h  cal¬ 
culated.  When  dealing  with  latitudes  greater  than  60° ,  however,  and  when  the  zenith 
angles  of  interest  are  90°  to  108°,  there  are  days  of  the  year  when  these  angles  do 
not  actually  occur;  that  is,  these  zenith  angles  are  not  reached.  In  the  generation  of 
the  tables,  these  angles  are  inserted  in  the  basic  formula  as  a  value  for  (z),  resulting 
in  a  computed  value  of  (h)  which  is  also  in  error.  Therefore,  the  time  computed  is 
also  in  error.  This  possibility  has  been  considered,  and  the  program  will  reject  all 
erroneous  times  that  are  calculated  due  to  the  insertion  of  an  angle  that  is  not  physi¬ 
cally  realizable;  and  insert  a  series  of  dash  marks  in  the  space  ordinarily  occupied 
by  the  time.  See  Appendix  D  for  a  sample  output  page  which  gives  twilight  informa¬ 
tion  for  Thule,  Greenland,  on  September  10,  1968. 


HI.  Use  of  the  Program 

This  section  is  intended  to  provide  all  the  information  necessary  to  apply  this 
program  to  any  first  generation  or  newer  computer  capable  of  being  programmed  in 
Fortran.  However,  the  input  and  output  statements  used  in  the  program  as  listed 
in  Appendix  B  may  have  to  be  modified  to  allow  its  use  on  any  specific  computer. 
The  body  of  the  program  may  be  used  as  listed. 

Input  data  and  format 

The  information  given  below  is  in  the  order  that  it  necessarily  must  be  placed 
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in  a  data  deck.  We  will  consider  as  an  example  the  generation  of  a  set  of  tables  for 
the  months  of  September  and  October,  1968,  at  a  location  of  58°  40'  latitude  and 
4$°  33'  longitude.  We  desire  the  times  associated  with  the  horizon  angles  from  -18° 
to  0°  in  1/2°  increments.  This  would  correspond  to  zenith  angles  of  108°  to  90°. 
Card  1.  Longitude  and  latitude,  with  west  longitude  and  north  latitude  being  taken 
as  positive. 

Enter  the  data  with  input  format  (2F9,  4).  Example, 

;  :  42.5500:  :  58.6667,  where  indicates  a  space. 

Card  2.  Increment  between  tabulated  angles;  maximum  zenith  angle,  minimum 
zenith  angle,  and  value  of  AT  as  taken  from  the  "AE  and  NA" 

Enter  the  data  with  the  input  format  (F4. 1,  2F5.  0,  F9.  4) 

Example,  :0.  5:1 08L:  :90.  ;37.  0000 
Card  3.  Month 

Enter  the  month  number,  1  thru  12,  with  the  input  format  (12). 

Example,  ;9. 

Cards  4+.  Day  of  month,  apparent  declination  in  degrees,  minutes,  and  seconds, 
and  ephemeris  transit  in  hours,  minutes,  and  seconds. 

Enter  the  data  with  the  input  format  (12,  15,  13,  F5. 1,  213,  F6.2). 
Example,  ;1  ;  ;  :18  ;  ;  3:53.4:12  ;  :  6:13.83. 
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0  22  50.7 

11 

52 

38.89 

23 

0-30-31.8 

11 

52 

18.00 

24 

0-23-55.2 

11 

51 

57.23 

25 

0-47-18.  9 

11 

51 

36.60 

26 

-1  10  42.7 

11 

51 

16.14 

34 

-4  177  9.  8 

11 

48 

40.48 

50 

10 

1 

-3  7  29.0 

11 

49 

36.  86 

35 

-15  1942.  '} 

11 

43 

36.  51 

99 

Notice  that  there  are  four  more  days  data  than  there  are  days  in  the  month. 

These  extra  days  are  the  first  four  days  of  the  following  month  and  are  necessary 
for  accurate  interpolation,  of  the  input  values  of  apparent  declination  and  ephemeris 
transit. 

Notice  also  the  manner  in  which  the  sign  of  the  apparent  declination  data 
changes.  Unless  the  degrees  term  is  an  integer  other  than  zero  when  the  minutes 
and  seconds  are  negative,  the  minutes  and  seconds  terms  must  be  given  the  minus  sign. 


3 


Control 

If  tables  are  to  be  generated  for  several  months  during  the  same  computer  run, 
the  last  entry  is  the  day  number  column  for  all  months  except  the  la^t  must  be  50, 
followed  by  the  next  monthfs  data.  This  allows  the  following  month's  data  to  be  fed 
in  automatically.  To  stop  the  program  at  the  end  of  the  last  month's  data,  enter  a 
99  in  the  day  number  column  using  the  same  format  as  the  day  number,  (12),  See 
Appendix  C  and  the  example  given  in  the  directions  for  cards  4  . 

It  should  be  obvious  that  there  is  no  need  to  generate  a  whole  month's  table  if  the 
information  is  needed  only  for  a  small  number  of  days.  For  example,  if  the  period 
of  interest  were  from  Sept.  8  to  Sept.  10,  1968,  it  would  be  necessary  only  to  insert 
ephemeris  data  for  Sept.  8  to  Sept.  14,  1968.  As  explained  previously,  the  extra 
four  days  data  are  used  for  the  purposes  of  interpolation. 
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Appendix  A 


Program  Flow  Chart 


start 


rread  long. ,  lat. ,  zenith  max  and  min  angles' 


read  month,  1st  five  day  ephemeris  data 


print  column  headings" 


compute  ephemeris  difference  table 


compute  interpolation  factor 


interpolate  ephemeris  data 


compute  hour  angle 


atime  :  i! 


atime  =  trise 


compute  trise 


trise  -  atime  :  .  ooj 


atime  =  trise  +12 


compute  GMT  of  trise  and  tset 


check  for  erroneous  times 


compute  tset 


jet  -  atime  :  .  00C 


print  mo. ,  day,  gmtrse,  gmtse 


^ang  ;  zrnim 


ang  =  ang  +  inc 


Appendix  B 
Program 

DIMENSION  IAD  1(5),  IAD2(5),  AD3(5),  APD(5),  AD)  (5),  AD2(5) 
DIMENSION  IER  1(5),  IER2(5),  ER3(5),  ETR(5),  ER1  (5),  ER2(5) 
DIMENSION  DAPD1  (4),  DAPD2(3),  DAPD3(2),  IDY(5) 

DIMENSION  DETR1  (4),  DETR2(3),  DETR3(2) 
c  THE  LONGITUDE  AND  LATITUDE  ARE  NOW  READ  IN. 

READ  7,  ALONG,  ALAT 
7  FORMAT  (2F9.  4) 

c  THE  ANGLE  INCREMENT,  ZENITH  MAXIMUM  AND  MINIMUM 

c  ANGLE  AND  VALUE  OF  DELTA  T  ARE  NOW  READ  IN. 

READ  9,  AINC,  ZANGMX,  ZANGMN,  DELTAT 

9  FORMAT  (F4.1,  2F5.  0,  F9.  4) 

c  DEGREE  TO  RADIAN  CONVERSION  FACTORS  AND  INITIAL 

c  VALUE  OF  APPROXIMATE  TIME  ARE  COMPUTED. 

ZANGMX  =  ZANGMX  +  AINC 
DTOR  =  3.  1415926/180. 

RTOD  =  180./3.  1415926 
ATIME  =  6.00 

c  THE  MONTH  NUMBER  AND  THE  FIRST  FIVE  DAYS 

c  DATA  ARE  READ  IN,  TRIG.  CONSTANTS  COMPUTED, 

c  AND  OUTPUT  HEADINGS  PRINTED. 

911  READ  912,  MO 

912  FORMAT(I2) 

DO  11  1=1,5 

II  READ  10,  IDY(I),  IAD1  (I),  IAD2(I),  AD3(I),  IER1  (I),  IER2(I),  FR3(I) 

10  FORMAT  (12,  15,  13,  F5.  1,  213,  F6.  2) 

ALATRD  =  ALAT*DTOR 

A  =  (SINF  (ALATRD))/  (COSF  (ALATRD)) 

E  =  1.  /COSF  (ALATRD) 

III  PRINT  890,  ALONG,  ALAT 

890  FORMAT  (1  HI,  9HLOCATION,  5X,  9HLONGITUDE,  FI  1.4,  5X, 

1  8HLATITUDE,  FI  1.4) 

PRINT  891 

891  FORMAT  (10X) 

PRINT  363 
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Program  (continued) 

363  FORMAT  (IX,  5HMONTH,  2X,  3HDAY,  4X,  7  HHORIZON,  5X, 

1  7HMORNING,  13X,  7HEVENING) 

PRINT  365 

365  FORMAT  (15X,  5HANGLE,  7X,  1 1HHR  MIN  SEC,  9X,  1 IHHR  MIN  SEC) 
PRINT  366 

366  FORMAT  (/) 

ANG  =  0.0 
SPACER  =  1.0 

c  ALL  DATA  IN  MINUTES  AND  SECONDS  ARE  CONVERTED 

c  TO  FRACTIONS  OF  HOURS  OR  DEGREES  AND  A  DIFFERENCE 

c  TABLE  IS  COMPUTED. 

12  DO  21  J  =  1,  5 

ADI  (J)  =  IAD1  (J) 

AD2  (J)  =  IAD2  (J) 

ER1  (J)  =  IER1  (J) 

ER2  (J)  =  IER2  (J) 

IF  (AD1(J)  )  18,  20,  20 

18  APD  (J)  =  ADI  (J)  -  (AD2  <J)/60.  )  -  (AD3  (J)/3600. ) 

GO  TO  21 

20  APD(J)  =  AD1(J)  +  (AD2(J)/60. )  +  (ADS', ,  ;/3600.  ) 

21  ETR(J)  =  ER1(J)  +  (ER2(J)/60. )  +  (ER3(J)/3600.  ) 

ro  26  K  =  1,4 

DAPD1  (K)  =  APD  (K+l)  -  APD  (K) 

26  DETR1  (K)  =  ETR  (K+l)  -ETR  (K) 

DO  30  L  =  1,  3 

DAPD2(L)  =  DAPDl(LH)  -  DAPDl(L) 

30  DETR2(L)  =  DETRKL+l)  -  DETRl(L) 

DO  34  M  =  1,2 

DAPD3(M)  =  DAPD2(M+1)  -  DAPD2(M) 

34  DETR3(M)  =  DETR2(M+1)  -  DETR2(M) 

DAPD4  =  DAPD3(2)  -  DAPD3(1) 

DETR4  =  DETR3(2)  -  DETR3(1) 

D  =  ZAMGMZ*DTOR 

c  AN  INTERPOLATION  FACTOR  IS  COMPUTED  AND  WITH 

c  IT  THE  APPARENT  DECLINATION  AND  EPHEMERIS  TRANSIT 

c  CORRESPONDING  TO  THE  APPROXIMATE  TIME  ARE  COMPUTED 

38  B  =  A  TIME/ 24.  +  1. 

B1  =  B-l. 

B2  =  B-2. 
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Program  (continued) 

B3  =  B-3. 

B4  =  B-4. 

ITLE  =  0 
ITLF  =  0 

AD  =  APD(l)  +  DAPD1(1)*B1  +  DAPD2(1  )*Bl*B2/2. 

1  +DAPD3(1  )*Bl*B2*B3/6.  +  DAPD4*Bl*B2*B3*B4/24. 

ET  =  ETR(l)  +  DETR1(1  )*B1  +  DETR2(l)*Bl*B2/2. 

1  +DETR3(l)*Bl*B2*B3/6.  +  DETR4*B1  B2*B3*B4/24. 

c  USING  THE  VALUES  OF  APPARENT  DECLINATION  AND  EPHEMERIS 

c  TRANSIT  COMPUTED  BY  INTERPOLATION,  THE  LOCAL  HOUR 

c  ANGLE  IS  COMPUTED  AND  CONVERTED  TO  AN  EQUIVALENT  TIME. 

ADRAD  =  AD*DTOR 
C  =  (SINF  (ADRA D))/ (COSF  (ADRAD)) 

FE  =  I.  /COSF  (ADRAD) 

G  =  COSF  (D) 

H  COS  =  -A*C+E*FE*G 
O  =  HCOS*  *2 
X  =  ABSF  (l.-O) 

HSIN  -  SQRTF  (X) 

ABHCOS  =  ABSF  (HCOS) 

HRAD  =  ATANF  (HTAN) 

IF  (HCOS)  59,  61,  61 

59  HDEGRE  =  180.  -  HRAD*RTOD 

GO  TO  62 

61  HDEGRE  =  KRAD*RTOD 

62  HTIME  -  HDEGRE/ 15. 

c  IT  IS  NEXT  DETERMINED  WHETHER  THE  APPROXIMATE  TIME 

c  USED  CORRESPONDS  TO  A  RISING  OR  SETTING  ANGLE  AND 

c  IF  THE  TOLERANCE  CONDITIONS  IS  MET.  IF  NOT,  THE 

c  APPROXIMATE  TIME  IS  REPLACED  BY  THE  COMPUTED  TIME  AND 

c  THE  ENTIRE  PROCESS  REPEATED  UN^IL  THE  TOLERANCE 

c  CONDITION  IS  MET. 

IF(12.  -ATIME)  70,  63,  63 

63  TRISE  =  ET-HTIME 

ABTIME  =  ABSFITRISE  -  ATIME) 

IF(.  0001  -  ABTIME)  66,  68,  68 
66  ATIME  =  TRISE 

GO  TO  38 

68  IF  (ANG)  661,  661,  663 
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Program  (continued) 


661  ATIME1  =  ATIME 

663  ATIME  =  TRISE  +  12. 

GO  TO  38 

70  TSET  =  ET  +  HTIME 

ACTIME  =  ABSF  (TSET  -  ATIME) 

IF  (.0001  -  ACTIME)  73,  75,  75 
73  ATIME  =  TSET 

GO  TO  38 

c  THE  COMPUTED  TIMES  ARE  CONVERTED  TO  GMT  AND  BROKEN 

c  DOWN  INTO  HOURS,  MINUTES  AND  SECONDS.  NEXT  A 

c  CHECK  IS  MADE  OF  THE  VALIDITY  OF  THE  TIMES  COMPUTED, 

c  AN  APPROPRIATE  OUTPUT  FORMAT  STATEMENT  SELECTED 

c  AND  THE  TIMES  PRINTED. 

75  GMTRSE  =  TRISE  =  ALONG/ 15.  -  DELTAT/3600. 

GMTSET  =  TSET  +  ALONG/ 15.  -  DELTAT/3600. 

IF(24.  -  GMTSET)  744,  745,  745 

744  GMTSET  =  GMTSET-24. 

745  ITRISE  =  GMTRSE 
ATRISE  =  ITRISE 

BTRISE  =  GMTRSE  -  ATRISE 
BTRISE  =  BTRISE*60. 

IBTRSE  =  BTRISE 

BTRSE  =  IBTRSE 

CTRISE  =  (BTRISE  -  BTRSE)*60. 

ITSET  =  GMTSET 
A  TSET  =  ITSET 
BTSET  =  GMTSET  -  A  TSET 
BSET  =  BTSET*60. 

IBSET  =  BSET 

A  BSET  =  IBSET 

CTSET  =  (BSET  -  ABSET)*60. 

IDAY  =  IDY(l) 

SPACER  =  SPACER  +  1.0 
IF  (2.0  -  SPACER)  779,  779,  800 

779  SPACER  =  0.0 
PRINT  780 

780  FORMAT  (10X) 

800  EDEP  =  (ZANGMX  -  ANG)  -90. 

ADEP  =  90. -(ZANGMX  -  ANG) 
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■I 


201 

202 

203 

204 

205 
806 

90 

807 

808 

809 

810 
811 
812 

813 

814 

815 


c 

c 

c 

c 

c 

c 

200 

91 

92 
930 


95 


Program  (continued) 

IF(ANG-(AINC/2. ))  200,  201,  201 
IF(TRISE  -  TRISE1 )  202,  203,  203 
ITLE  =  5 

IF(TSET1  -  TSET)  204,  205,  205 
ITLF  =  10 

IF(ITLE  +  ITLF)806,  806,  807 
PRINT  90,  MO,  IDAY,  ADEP,  ITRISE, 

IBTRSE,  CTRISE,  ITSET,  IBSET,  CTSET 

FORMAT  (IX,  2(2X,  13)  F8.  1,  5X  214,  F6.  1,  6X,  214,  F6.1) 

GO  TO  200 

IF  ((ITLE  +  ITLF)  -5)  200,  808,  810 

PRINT  809,  MO,  IDAY,  ADEP,  ITSET,  IBSET,  CTSET 

FORMAT  (IX,  2(2X,  13),  F8.  1,  7X,  12H . ,214,  F6.1) 

GO  TO  200 

IF  ((ITLE  +  ITLF)-10)  200,  811,  813 

PRINT  812,  MO,  IDAY,  ADEP,  ITRISE,  IBTRSE,  CTRISE 

FORMAT  (IX,  2(2X,  13).  F8.  1,  5X,  214,  F6.1,  8X,  12H--  --  - ) 

GOTO  200 

IF  ((ITLE  +  ITLF)- 15)  200,  814,  814 
PRINT  815,  MO,  IDAY,  ADEP 

FORMAT  (IX,  2(2X,  13),  F8.  1,  7X,  12H . -, 

8X,  12H--  --  - ) 

NEXT  IT  IS  DETERMINED  IF  ALL  THE  ANGLES  FOR  THE  DAY 
HAVE  BEEN  COMPUTED.  IF  NOT,  THE  PROGRAM  INCREMENTS 
TO  THE  NEXT  ONE  AND  THE  TIMES  ARE  COMPUTED. 

IF  ALL  THE  ANGLES  FOR  THE  DAY  HAVE  BEEN  COMPUTED, 

THE  PROGRAM  INCREMENTS  TO  THE  NEXT  DAYS  DATA 
AND  CONTINUES  AS  BEFORE. 

TRISE1  =  TRISE 
TSET1  =  TSET 

IF  (ZANGMN  -  (EDEP  +  90.  ))92,  95,  95 

D  =  D-AINC*DTOR 

ANG  =  ANG  +  AINC 

A  TIME  =  TRISE 

GO  TO  38 

ATIME  =  TRISE 

DO  102  J  1,4 

IDY(J)  =  IDY(J+1 ) 

IADl(J)  =  IAD1  (J+l) 
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Program  (continued) 


IAD2(J)  =  IAD2(J+1 ) 

AD3(J)  =  AD3(J+1 ) 

IERl(J)  =  IERKJ+1) 

IER2(J)  =  IER2(J+l ) 

102  ER3(J)  =  ER3(J+1) 

READ  10,  IDY(5),  IAD1(5),  IAD2(5),  AD3(5  ),  IER1(5),  IER2(5),  ER3(5) 
IF  (IDY(5)  -  50)  104,  125,  126 
104  SPACER  =  -1.0 

GO  TO  111 

125  GO  TO  911 

126  STOP 
END 
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Sample  Ephemeris  Input  Data  (continued) 
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1 

-3 

7 

29. 

0 

11 

49 

36. 

86 

2 

-3 

30 

45. 

2 

11 

49 

17. 

76 

3 

-3 

53 

58. 

9 

11 

48 

58. 

96 

4 

-4 

17 

9. 

8 

11 

48 

40. 

48 

5 

-4 

40 

17. 

5 

11 

48 

22. 

35 

6 

-5 

3 

21. 

8 

11 

48 

4. 

60 

7 

-5 

26 

22. 

2 

11 

47 

47. 

24 

8 

-5 

49 

18. 

6 

11 

47 

30. 

29 

9 

-6 

12 

10. 

4 

11 

47 

13. 

79 

10 

-6 

34 

57. 

5 

11 

46 

57. 

75 

11 

-6 

57 

39. 

4 

11 

46 

42. 

20 

12 

-7 

20 

15. 

8 

11 

46 

27. 

14 

13 

-7 

42 

46. 

3 

11 

46 

12. 

62 

14 

-8 

5 

10. 

5 

11 

45 

58. 

63 

15 

-8 

27 

28. 

1 

11 

45 

45. 

21 

16 

-8 

49 

38. 

7 

11 

45 

CO 

to 

37 

17 

-9 

11 

41. 

9 

11 

45 

20. 

13 

18 

-9 

33 

37. 

2 

11 

45 

9. 

50 

19 

-9 

55 

24. 

4 

11 

44 

57. 

50 

20 

-10 

17 

2. 

9 

11 

44 

47. 

14 

21 

-10 

38 

32. 

4 

U 

44 

37. 

45 

22 

-10 

59 

52. 

4 

11 

44 

28. 

41 

23 

-11 

21 

2. 

6 

11 

44 

20. 

06 

24 

-11 

42 

2. 

4 

J1 

44 

12. 

39 

25 

-12 

2 

51. 

6 

11 

44 

5. 

43 

26 

-12 

23 

29. 

6 

11 

43 

59. 

17 

27 

-12 

43 

56. 

1 

11 

43 

53. 

62 

28 

-13 

4 

10. 

5 

11 

43 

48. 

81 

29 

-13 

24 

12. 

6 

11 

43 

44. 

73 

30 

-13 

44 

1. 

8 

11 

43 

41. 

41 

31 

-14 

3 

37. 

9 

11 

43 

38. 

85 

32 

-14 

23 

0. 

3 

11 

43 

37. 

07 

33 

-14 

42 

8. 

8 

11 

43 

36. 

08 

34 

-15 

1 

3. 

0 

11 

43 

35. 

89 

35 

-15 

19 

42. 

5 

11 

43 

36. 

51 

99 


I 


! 

I 
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Location, 

Month 


Appendix  D 


Sample  Output 


Longitude  68. 6667 


76.5500 


Day 


Horizon  Morning 

Angle  Hr  Min  Sec 


Evening 
Hr  Min  Sec 


Latitude 


* 


Appendix  E 
Sample  Calculation 

Reproduction  of  example  from  section  13.  D  of  the  Explanatory  Supplement 
to  the  American  Ephemeris  and  Nautical  Almanac. 


March  7,  1960 

Latitude  52° 

i  =  52° 

tan 

i  =  1.2799 

sec 

i  =  1.6243 

Morning 

Twilight 

Sunrise 

Sunset 

Evening 

Twilight 

Approx,  time 

4h  30m 

6h  30m 

1 8h  00m 

19h  30m 

5  (A.  E. ,  pg  21) 

-5°  17' 

-5°  15' 

-5°  04' 

-5°  02' 

tangent  6 

-0. 0925 

-0.  0919 

-0. 0887 

-0.0881 

secant  6 

1. 0043 

1. 0042 

1. 0039 

1. 0039 

z 

108° 

CD 

O 

O 

C Jl 

o 

90°  50' 

108° 

cosine  z 

-0.  3000 

-0.  0145 

-0. 0145 

-0.3090 

-tan  $tan  6 

0.  1184 

0.  1176 

0.  1135 

0. 1128 

cos  z  sec  $  sec  6 

-0. 5041 

-0.0237 

-0. 0237 

-0.5039 

cosine  h 

-0. 3857 

-0. 0939 

0. 0898 

-0.3911 

h  min 

h  min 

h  min 

h  min 

h 

7  30.7 

5  38.5 

5  39.4 

7  32.1 

ephemeris  transit 

12  11.? 

12  11.2 

12  11.0 

12  11.0 

Time  (local) 

4  40.  5 

6  32.7 

17  50.4 

19  43.1 
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Appendix  F 


IAD1 

IAD2 

AD3 

A  PD 

ADI 

AD2 

IER1 

IER2 

ER3 

ETR 

ER1 

ER2 

DAPD1 

DAPD2 

DAPD3 

DAPD4 

DETR1 

DETR2 

DETR3 

DETR4 

IDY 

ALONG 

ALAT 

AINC 

ZANGMX 

ZANGMN 

DELTA T 

DTOR 

RTOD 

ATIME 

MO 

ALATRD 

A 

E 

ANG 


Glossary  of  Terms 

apparent  declination,  degrees  (fixed  point) 
apparent  declination,  minutes  (fixed  point) 
apparent  declination,  seconds 

apparent  declination,  degrees  and  decimals  of  degrees 
apparent  declination  ,  degrees  (floating  point) 
apparent  declination,  minutes  (floation  point) 
ephemeris  transit,  hours  (fixed  point) 
ephemeris  transit,  minutes  (fixed  point) 
ephemeris  transit,  seconds 

ephemeris  transit,  hours  and  decimals  of  hours 
ephemeris  transit,  hours  (floating  point) 
ephemeris  transit,  minutes  (floating  point) 
apparent  declination,  first  difference 
apparent  declination,  second  difference 
apparent  declination,  third  difference 
apparent  declination,  fourth  difference 
ephemeris  transit,  first  difference 
ephemeris  transit,  second  difference 
ephemeris  transit,  third  difference 
ephemeris  transit,  fourth  difference 
day  number 
longitude 
latitude 

degree  increment 
maximum  zenith  angle 
minimum  zenith  angle 

correction  factor,  ephemeris  time  to  universal  time 
conversion  factor,  degrees  to  radians 
conversion  factor,  radians  to  degrees 
approximate  time 
month 

latitude  in  radians 
tangent  (latitude) 
secant  (latitude) 
counter  for  angle  increment 
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Glossary  of  Terms  (continued) 


SPACER 

D 

B 

B1  TO  B4 

ITLE 

ITLF 

AD 

ET 

ADRAD 

C 

FE 

G 

HCOS 

HSIN 

ABHCOS 

HTAN 

HRAD 

HDEGRE 

HTIME 

TRISE 

ABTIME 

TSET 

ACTIME 

GMTRSE 

GMTSET 

ITRISE 

ATRISE 

BTRISE 

IBTRSE 

BTRSE 

CTRISE 

ITSET 

ATSET 

BTSET 

IBSET 

BSET 

ABSET 

CTSET 


counter  to  control  output  spacing 

zenith  angle  in  radians 

interpolation  factor 

interpolation  factors 

error  check  counter 

error  check  counter 

apparent  declination  after  interpolation 

ephemeris  transit  after  interpolation 

apparent  declination  in  radians 

tangent  (apparent  declination) 

secant  (apparent  declination) 

cosine  (zenith  angle) 

cosine  of  hour  angle 

sine  of  hour  angle 

absolute  value  of  the  cosine  of  the  hour  angle 
tangent  of  hour  angle 
hour  angle  in  radians 
hour  angle  in  degrees 
time  corresponding  to  hour  angle  of  H 
local  time  of  zenith  angle  during  AM 
absolute  value  of  local  time  during  AM 
local  time  of  zenith  angle  during  PM 
absolute  value  of  local  time  during  PM 
GMT  of  hour  angle  during  AM 
GMT  of  hour  angle  during  PM 
hour  of  AM  zenith  angle  (fixed  point) 
hour  of  AM  zenith  angle  (floating  point) 
minutes  of  AM  zenith  angle  (floating  point) 
minutes  of  AM  zenith  angle  (fixed  point) 
minutes  of  AM  zenith  angle  (floating  point) 
seconds  of  AM  zenith  angle  (fixed  point) 
hours  of  PM  zenith  angle  (fixed  point) 
hours  of  PM  zenith  angle  (floating  point) 
minutes  of  PM  zenith  angle  (floating  point) 
minutes  of  PM  zenith  angle  (fixed  point) 
minutes  of  PM  zenith  angle  (floating  point) 
minutes  of  PM  zenith  angle  (floating  point) 
seconds  of  PM  zenith  angle  (floating  point) 


Glossary  of  Terms  (continued) 


IDAY 

day  for  rrintout 

EDEP 

counter  for  zenith  angle 

ADEP 

counter  for  zenith  angle 

TRISE1 

ne  of  zenith  angle  AM 

TSET1 

time  of  zenith  angle  PM 
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